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ABSTRACT 

Turbulence dynamo deals with amplification of a seed magnetic field in a turbulent medium and 
has been studied mostly for uniform or spatially homogeneous seed magnetic fields. However, some 
astrophysical processes (e.g. jets from active galaxies, galactic winds, or ram-pressure stripping in 
galaxy clusters) can provide localized seed magnetic fields. In this paper, we numerically study 
amplification of localized seed magnetic fields in a turbulent medium. Throughout the paper, we 
assume that driving scale of turbulence is comparable to the size of the system. Our findings are as 
follows. First, turbulence can amplify a localized seed magnetic field very efficiently. The growth rate 
of magnetic energy density is as high as that for a uniform seed magnetic field. This result implies 
that a magnetic field ejected from an astrophysical object can be a viable source of magnetic field in 
a cluster. Second, the localized seed magnetic field disperses and fills the whole system very fast. If 
turbulence in a system (e.g. a galaxy cluster or a filament) is driven at large scales, we expect that it 
takes a few large-eddy turnover times for magnetic field to fill the whole system. Third, growth and 
turbulence diffusion of a localized seed magnetic field are also fast in high magnetic Prandtl number 
turbulence. Fourth, even in decaying turbulence, a localized seed magnetic field can ultimately fill the 
whole system. Although the dispersal rate of magnetic field is not fast in purely decaying turbulence, 
it can be enhanced by an additional forcing. 

Subject headings: ISM:general - intergalactic medium - MHD - turbulence 



1. INTRODUCTION 

The origin and growth of magnetic field in the universe 
is an important unsolved problem in astrophysics (Kul- 
srud & Zweibel 2008). Earlier studies have shown that 
turbulence can efficiently amplify a weak seed magnetic 
field. The main mechanism that is responsible for the 
amplification is stretching of magnetic field lines (Batch- 
clor 1950; Zel'dovich et al. 1984; Childress & Gilbert 
1995). This process is called turbulence dynamo (or 
small-scale dynamo) and should not be confused with 
mean field dynamo (or large-scale dynamo), which deals 
with generation/growth of a large-scale magnetic field. 

Small-scale turbulence dynamo has been studied since 
1950s (see, e.g., Batchelor 1950; Kazantsev 1968; Kul- 
srud & Anderson 1992; Kulsrud et al. 1997; Cho & Vish- 
niac 2000; Schekochihin et al. 2004; Brandenburg & Sub- 
ramanian 2005; Schekochihin & Cowley 2007; Ryu et al. 
2008; Cho et al. 2009). When we introduce a weak 
uniform magnetic field in a turbulent medium, amplifi- 
cation of the uniform field happens in three stages (see 
Schekochihin & Cowley 2007; see also Cho & Vishniac 
2000, Cho et al. 2009). (1) Stretching of magnetic field 
lines occurs most actively near the velocity dissipation 
scale first, and the magnetic energy grows exponentially. 
Note that eddy motions in turbulence are fastest at the 
velocity dissipation scale. (2) The exponential growth 
stage ends when the magnetic energy becomes compara- 
ble to the kinetic energy at the dissipation scale. The 
subsequent stage is characterized by a linear growth of 
magnetic energy and a gradual shift of the peak of mag- 
netic field spectrum to larger scales. (3) The amplifi- 
cation of magnetic field stops when the magnetic energy 
density becomes comparable to the kinetic energy density 
and a final, statistically steady, saturation stage begins. 



Knowledge on turbulence dynamo is useful for under- 
standing the origin of magnetic fields in many astrophys- 
ical fluids. For example, Ryu et al. (2008) numerically 
studied amplification of the intergalactic magnetic field 
by the dynamo action of turbulence. They assumed that 
turbulence is generated by cosmological shocks, which 
are produced when primordial gases fall into the grav- 
itational potential well caused by dark matters in clus- 
ters and filaments (Ryu et al. 2003). For simplicity, 
they assumed that a weak uniform seed magnetic field 
fills the whole system at the beginning of the simulation. 
They found that turbulence is strong in the intraclus- 
tcr medium (ICM) and, therefore, has almost reached 
the saturation stage and that turbulence in filaments is 
weak and in early linear stage. Based on a more quanti- 
tative analysis, they estimated that the average magnetic 
field strength would be a few /iG inside clusters/groups, 
approximately 0.1 /iG around clusters/groups, and ap- 
proximately 10 nG in filaments. 

Existence of regular magnetic fields in many astrophys- 
ical fluids, such as the interstellar medium (ISM), justi- 
fies the use of uniform seed magnetic fields in turbulence 
dynamo studies. However, there are astrophysical sys- 
tems where existence of regular fields is uncertain. For 
example, it is not clear whether or not regular fields 
exist in the ICM. If seed magnetic fields are relic pri- 
mordial fields (Rees 1987; Kronberg 1994; Dolag, Bartel- 
mann & Lesch 1999, 2002; Gnedin, Ferrara & Zweibel 
2000; Banerjee & Jedamzik 2003; Murgia et al. 2004; 
Briiggen et al. 2005; Dubois & Teyssier 2008), the cor- 
relation length of the seed fields can be very large and 
we may have regular fields in the ICM. On the other 
hand, if seed magnetic fields have astrophysical origins 
(e.g. jets from active galaxies, galaxy winds, and ram- 
pressure stripping, etc.; see Rephaeli 1988 for an ear- 
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lier discussion on the origin of the ICM magnetic fields 
through ram-pressure stripping), we will probably have 
no regular fields in the ICM. In the latter case, we ex- 
pect that seed magnetic fields ejected from astrophysical 
objects are highly localized in space. Then a question 
arises: how do these localized seed magnetic fields grow? 
In this paper, we study amplification of localized seed 
magnetic fields in a turbulent medium. 

Earlier cosmological simulations have shown that mag- 
netic fields ejected from active galactic nuclei (AGNs) 
(Ensslin et al. 1997; Volk & Atoyan 2000; Donnert et al 
2009; Xu et al. 2010, 2011) or normal galaxies through 
galactic winds (Kronberg, Lesch & Hopp 1999; Bertone, 
Vogt & Ensslin 2006) or ram-pressure stripping (Arieli, 
Rephaeli & Norman 2011) can be dispersed and amplified 
by turbulence and are viable origins of magnetic fields in 
galaxy clusters. Xu et al. (2010), for example, showed 
that "as long as the AGN magnetic fields are ejected be- 
fore the major mergers in the cluster formation history, 
magnetic fields can be transported throughout the clus- 
ter and can be further amplified by the ICM turbulence 
caused by hierarchical mergers during the cluster forma- 
tion process." However, those simulations are compress- 
ible ones and it is usually difficult to obtain converging 
results free from numerical effects in compressible tur- 
bulence simulations (see discussions in Ryu et al. 2008). 
Therefore it is advantageous to use an incompressible 
code because we can control numerical viscosity and re- 
sistivity in incompressible codes. 

In this paper, we use an incompressible code with 
hyper-viscosity and hyper-diffusion, so that we have vir- 
tually zero intrinsic numerical viscosity and diffusion, 
and the maximized inertial range. We assume that driv- 
ing scale of turbulence is comparable to the size of the 
whole system, which would be true if turbulence is pro- 
duced by cosmological shocks or major mergers in galaxy 
clusters. We also assume that the seed magnetic field has 
a doughnut shape at t=0. We describe numerical method 
in Section 2, and we present results in Section 3. We give 
discussions in Section 4 and the summary in Section 5. 

2. NUMERICAL METHODS 

We use a pseudospectral code to solve the incompress- 
ible MHD equations in a periodic box of size 27r: 

Ov 

— = -(Vxv)xv+(VxB)xB + iA7 2 v + f+VP', (1) 
ot 

— =Vx (vxBji/B, (2) 

V • v = V • B = 0, (3) 

where f is random driving force, P' = P + v- v/2, v is 
the velocity, and B is magnetic field divided by (4-Kp) 1 / 2 . 
We use 22 forcing components with 2 < k < \/l2. Each 
forcing component has correlation time of ~ 1. The peak 
of energy injection occurs at k rs 2.5. In our simulations 
with v — r\, v ~ 1 during the growth stage and ~ 0.8 
during the saturation stage (see Figure 1). Therefore, 
one eddy turnover time at the outer scale of turbulence 
(or large-eddy turnover time), ~ L/v, is approximately 
^2.5 and ^3 time units, respectively. Since v ~ 1, v 
can be viewed as roughly the velocity measured in units 
of the r.m.s. velocity of the system and B as the Alfven 



speed in the same units. Other variables have their usual 
meaning. 

In pseudospectral methods, we calculate the temporal 
evolution of the equations ([T| and ^ in Fourier space. 
To obtain the Fourier components of nonlinear terms, 
we first calculate them in real space, and transform back 
into Fourier space. We use exactly same forcing terms 
for all simulations. We use an appropriate projection 
operator to calculate VP' term in Fourier space and also 
to enforce divergence-free condition (V • v = V • B = 0). 
We use up to 512 3 collocation points. 

At t = 0, the magnetic field has a doughnut shape, 
which mimics a magnetic field ejected by a galactic out- 
flow^ We use the following expression for the magnetic 
fieldat t=0: 

B(r±,Ax) = -%^rie-^/ 2CT oV A * 2 / 8 ^, (4) 
2cr Q e 

where er = 4y/2, r± = {Ay 2 + Az 2 ) 1 / 2 , and Ax, Ay, 
and Az are distances measured from the center of the 
numerical box in grid units. The unit vector 9± is per- 
pendicular to (Ax, 0,0) and (0, Ay, Az). Note that the 
maximum strength of the magnetic field at t=0 is B max . 
Since cto = 4\/2 in Equation (UJ), the size of the magne- 
tized region at t=0 is ~ 16 in grid units, which is ~ 1/16 
of the whole system for runs with 256 3 grid points and 
~ 1/32 for 512 3 . Therefore, in a cluster of size^ IMpc, 
the size of the initially magnetized region corresponds to 
~ GOkpc and ~ 30fcpc, respectively. 

Hyperviscosity and hyperdiffusion are used for dissi- 
pation terms. The power of hyperviscosity is set to 3, 
such that the dissipation term in the above equation is 
replaced with 

-^(V 2 )"v, (5) 

where n = 3. The same expression is used for the mag- 
netic dissipation term. We list parameters used for the 
simulations in Table [I] Diagnostics of our code can be 
found in Cho & Vishniac (2000). 

3. RESULTS 

3.1. Turbulence with unit magnetic Prandtl number 
(v = r,) 

Figure [l] shows time evolution of v 2 and B 2 for different 
values oiB max . The initial magnetic field has a dough- 
nut shape (see Equation [4]) and B max is the maximum 
strength of the initial magnetic field. As in the uniform 
seed magnetic field cases, growth of magnetic energy den- 
sity occurs in three stages. Magnetic energy grows ex- 
ponentially first. We can see the exponential growth in 
the lower panel of Figure [II where we use a logarithmic 
scale for the y-axis. WhenB max is small (see, for exam- 
ple, 256--B maa; 0.001 and 256--B maa; 0.01), the growth rate 
is virtually independent of the value of B max . Then a 
slower linear growth stage follows, which we can clearly 
observe in the upper panel of the figure. All runs show 
similar growth rates in this stage. Finally, after t ~ 20 
in units of L/v (~ 2.5 in code time units), a statistically 

1 Note however that we use a weak uniform seed magnetic 
field for the reference runs REF256 and REF256-Pr, in which the 
strength of the mean field, Bo, is 0.0001 and we have only the mean 
field at t=0. 
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Time [ = t code /(L/v)] Time [ = t code /(L/v)] 

Fig. 1. — Time evolution of B 2 and v 2 in runs with v = rj. Each run starts from a fully developed turbulent velocity field and a localized 
seed magnetic field. The localized seed magnetic fields have a doughnut shape. All runs have the same numerical resolution. Different 
lines denote different values of B max . Time is given in units of the large-eddy turnover time defined by L/v: t = t co de/ (L/v), where t co[ j e 
is time in code units, L (~ 2.5) is the outer-scale of turbulence and v (~ 1) is the r.m.s. velocity before the saturation stage. Results arc 
from 256- Bmax O.OOf , 256-B ma x0.01, 256-B ma xO-l, and 2§6-B max l. Note the logarithmic scale for the vertical axis in the lower panel. 

Fig. 2. — Resolution study. We compare results of 256-B ma a;0.0f and 5f 2-S maI 0.0f . Two runs have different numerical resolutions, but 
other numerical setups are the same. Two results show good agreements. Note the logarithmic scale for the vertical axis in the lower panel. 




Time Hcode/(L/v)] k 

Fig. 3. — Comparison with the case of a uniform seed magnetic field. We compare the results of 256--B ma i0.01 and REF256. Note that 
we use a localized seed magnetic field in the former and a uniform seed field (Bq = 0.0001) in the latter. We find a good correspondence 
between two results. 

Fig. 4. — Time evolution of magnetic spectrum. The initial seed magnetic field (see the solid curve) looses its original shape quickly and 
forms a spectrum that has a peak near the dissipation scale (see the dotted curve). Then the spectrum goes up without changing the shape 
much (see the dashed line). At this point, the seed magnetic field already fills the most of the simulation box (see Figure [5] and Figure 
|SJ. After magnetic field fills the whole volume, time evolution of magnetic spectrum and energy density should be very similar to those of 
weak uniform seed magnetic field cases. Time is given in units of L/v. From Run 512-f? ma:E 0.01. 
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stationary saturation stage is reached. The saturation 
level of magnetic energy density is comparable to the ki- 
netic energy density. Figure [T] shows that the stronger 
the seed magnetic field (or B max ) is, the earlier the linear 
growth stage begins and the sooner the saturation stage 
is reached. 

Resolution study shows that the linear growth rate 
does not strongly depend on the numerical resolution. 
In Figure |2j we compare Runs 512-B ma:E 0.01 and 256- 
B max 0.01, in which the shape of seed magnetic field (see 
Equation [4] ) and numerical setup are identical and only 
the numerical resolutions are different. As we can see 
in the figure, the overall results are very similar. That 
is, the growth rates at the linear growth stage and the 
saturation levels of the magnetic energy density are very 
similar. In the lower panel of the figure, we plot the 
same quantities, but we use a logarithmic scale for the y- 
axis. The lower panel shows the growth rate of Run 512- 
BmaicO-01 during the exponential growth stage is higher 
than that of 256-B maa; 0.01, which is due to the fact that 
the eddy turnover time at the velocity dissipation scale 
of Run 512-B maa; 0.01 is shorter. 

According to Figure [T] the growth rate of magnetic 
energy during the linear growth stage seems to be inde- 
pendent of the value of B max . The linear growth rate in 
Figure [T] is also in good correspondence with that of uni- 
form seed magnetic field cases. When we compare the 
linear grow rates of Run 256-B max 0.01 (localized seed 
magnetic field) and Run REF256 (uniform seed magnetic 
field), they are very similar (Figure [3]). The saturation 
levels of the magnetic energy density m both simulations 
are also very similar. 

Three stages of magnetic energy growth are also sup- 
ported by the time evolution of magnetic spectrum. Fig- 
ure [4] shows how magnetic spectrum Eb(k) of Run 512- 
B max 0.01 evolves with time. When t < 10, stretching 
happens near the dissipation scale and, therefore, mag- 
netic energy spectrum peaks there. As magnetic field 
is amplified by stretching, the magnetic spectrum moves 
upward without changing its shape much (compare, for 
example, spectra for t=0.24 and t=2.4). When mag- 
netic spectrum 'touches' the velocity spectrum E v (k) at 
the dissipation scale, the exponential growth ends and 
the slower linear growth stage follows. Note that, when 
magnetic spectrum 'touches' the velocity spectrum at the 
dissipation scale, the magnetic energy density is approx- 
imately ~ (fcrf/fci) 2 ^ 3 times the kinetic one, where kd is 
the dissipation scale wavenumber and kh ~ 2.5, and it 
becomes 'visible' in Figure[2j The long dashed curves are 
the spectra of velocity and magnetic field during the sat- 
uration stage. Velocity spectrum E v (k) peaks at the en- 
ergy injection scale and magnetic spectrum Et,(k) peaks 
at a smaller scale. The logarithmic slope of the veloc- 
ity spectrum during the saturation stage is steeper than 
—5/3 for k < 10 and gets shallower for k > 10. The 
behavior of magnetic energy spectrum is very similar to 
that of a uniform seed field case (see, for example, Cho 
et al. 2009). 

So far, we have found surprising similarities between 
localized and uniform seed magnetic fields. They both 
follow similar growth stages, show similar growth rates at 
the linear growth stage, reach similar levels of magnetic 
energy saturation, and exhibit similar time evolution of 



magnetic spectrum. Then why are they all similar? The 
answer is fast magnetic diffusion. Turbulence enhances 
diffusion processes. Figure [5] visualizes fast magnetic dif- 
fusion in Run 256-B maa: 0.01. The localized initial mag- 
netic field (left panel) spreads out by turbulent motions. 
Note that it takes only ~ 2.4 large-eddy turnover times 
for the localized initial magnetic field to fill virtually the 
whole box (see the right panel) . The large-eddy turnover 
time is defined as L/v, where L (~ 2.5) is the outer scale 
and v (~ 1) the r.m.s. velocity. Therefore, turbulence 
makes the localized seed magnetic field fill the whole sys- 
tem within less than ^3 large-eddy turnover times. After 
magnetic field fills the whole box, time evolution of mag- 
netic energy density should be very similar to that of 
uniform or spatially homogeneous magnetic field cases 
(see, for example, Cho et al. 2009). 

Figure [6] confirms fast diffusion of magnetic field, in 
which we plot the standard deviation, a, of magnetic 
field distribution: 



r 2sl/2 



\B(x,t)\ 2 d 3 x 



f \B(x,t)\ 2 d 3 x 
J Xi\B(x,t)\ 2 d 3 x 



J \B(x,t)\ 2 d 3 x 



(6) 
(7) 

(8) 



where i =x, y, and z and a shown in the plot is normal- 
ized by the box size. Initially the standard deviation rises 
very quickly, which is due to diffusion of magnetic field 
within an energy-containing eddy. The duration of the 
quick rise is about 1.5 large-eddy turnover times, if we 
define the large-eddy turnover time as L/v. The slower 
growth after t/(L/v) ~ 1.5 is due to a slower diffusion 
of magnetic field on scales larger than the outer scale 
of turbulence. The behavior of a is not sensitive to the 
value of B max . 

Richardson's law states that the mean square separa- 
tion of two passive particles, a m , satisfies 



(9) 



da m 
~dt 

Therefore, we may write 

a{t) 2 / 3 - a(t ) 2 / 3 cx t - t Q . (10) 

(see, for example, Lesieur 1990; see also Cho et al. 
2003). The inset shows a 2 / 3 of magnetic field distribution 
roughly follows the Richardson's law when t/(L/v) < 1.5. 

3.2. Turbulence with a high magnetic Prandtl number 
(v » 77J 

The viscosity in galaxy clusters may not be negligi- 
ble (Schekochihin et al. 2004; Ruszkowski et al. 2004; 
Reynolds et al. 2005; Subramanian, Shukurov & Haugen 
2006) , while the magnetic diffusivity is still very small. In 
fact, if we use the Spitzer (1962) formula for the viscosity, 
the Reynolds number Lv/v, where L is the outer scale 
of turbulence in a cluster, v the r.m.s. velocity, and v 
the viscosity, is less than ~ 10 3 . Therefore, the magnetic 
Prandtl number (y/r\) may be very large in the ICM. 
In this subsection, we consider growth of localized seed 
magnetic fields in turbulence with non-negligible viscos- 
ity. The Reynolds number (Lv/v) we use in this section 
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Fig. 5. — Turbulence diffusion of magnetic field. From left to right panels, snapshots at t=0, 1.2, and 2.4 are shown. The seed magnetic 
field has a doughnut shape at t=0. At t = 2.4, the localized seed field fills the most of the computational box. Here time (t) is measured 
in units of L/v. From 256--B ma: c0.01. 
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Fig. 6. — Time evolution of the standard deviation of magnetic 
field distribution. The standard deviation is normalized by the box 
size, so that it becomes 0.5 when a homogeneous magnetic field fills 
the whole box. The standard deviation rises quickly when t < 1.5 
large-eddy turnover times, which happens when the localized mag- 
netic field fills the whole energy-containing eddy. After t ~ 1.4 
large-eddy turnover times, the diffusion rate of the magnetic field 
slows down. It takes ~ 3 large-eddy turnover times for the mag- 
netic field to fill the whole box. Results are from 256-B maa! 0.001, 
256-B ma!c 0.01, 256-B ma!C l, and 512-B m<M 0.01. 

is 100. The magnetic diffusivity is negligibly small 
because we use a hyper-diffusion. 

Figure J7|a) shows time evolution of B 2 and v 2 . Due to 
non-negligible viscosity in those simulations, small-scale 
velocity is strongly damped and the inertial range of tur- 
bulence is poorly resolved. Indeed the velocity spectrum 
at t=0.0693 (solid curve) in Figure [Wb) confirms this. 
Therefore, stretching of magnetic field lines happens at 
a scale very close to the outer scale of turbulence and 
the growth rate of magnetic energy density is mostly ex- 
ponential (see Figure fflc); see also next paragraph for 
further discussions). Figure [7^ a) shows that it takes less 
than ~15 large-eddy turnover times for the runs to reach 
the saturation stage. Comparing with Figure [T] we note 
that growth of magnetic field is at least as fast as that 
of unit magnetic Prandtl number cases. Note that the 
large-eddy turnover time L/v in high magnetic Prandtl 



number cases is approximately 2.6 in code units before 
saturation and 4.5 in code units after saturation. 

It seems that growth of a localized seed magnetic field 
also happens in 3 stages in a high magnetic Prandtl num- 
ber turbulenc^Hl 

1. First, we can clearly observe the exponential 
growth stage in Figure ff\c). The exponential 
growth stage ends when the magnetic energy den- 
sity becomes comparable to the kinetic energy den- 
sity at the stretching scale. Since the stretching 
scale is very close to the outer scale of turbulence 
in our simulations, the exponential stage ends when 
magnetic energy density becomes comparable to 
the total kinetic energy density. Therefore, the ex- 
ponential stage takes most of the growth time in 
the high magnetic Prandtl number turbulence. Ve- 
locity spectrum does not change much during the 
exponential stage (compare E v (k) , s at t=0.693 and 
8.54 in Figure mb)), which means that magnetic 
back-reaction is negligible during the exponential 
growth stage. 

2. Second, after the exponential growth stage, a 
slower growth stage follows. In Run 512-B ma:r 0.01- 
Pr, the slower growth stage is relatively short due 
to proximity between the outer scale and the ini- 
tial stretching scale: it begins at t~12 and ends at 
t~17. According to Figure (7^a), the growth rate 
during the slower growth stage is roughly linear. 
The growth rate of B 2 during the slower growth 
stage can be estimated as follows. If stretching 
happens at a scale /, we can write 



(see Schekochihin & Cowley 2007), where vi is the 
r.m.s. velocity at the scale I. The stretching scale I 
increases during the slower growth stage. The peak 
of the magnetic spectrum in Figure \j\b) moves to 
smaller wavenumbers (compare magnetic spectra 
at t=12.69, 13.62, and 20.08), which might indi- 

2 For a related discussion on the case of a weak uniform seed 
magnetic field, see, for example, Schekochihin & Cowley (2007). 
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cate increase of the stretching scale. The veloc- 
ity spectrum between k ^3 and k ~10 becomes 
quickly suppressed as soon as the slower growth 
stage begins (compare velocity spectra at t=12.69 
and t=13.62). The wavenumbers at which velocity 
spectrum is suppressed might have something to 
do with the stretching scale during the exponential 
growth stage. After the suppression, the stretching 
scale is clearly the outer scale of turbulence (see the 
shape of J5„(fc)'s at t=13.62 and 20.08). Therefore, 
the behavior of velocity spectrum also suggests in- 
crease of I. We also note that decrease of v is not 
negligible during the slower growth stage: v 2 drops 
fast after t ~12 in Run 512-B maa: 0.01-Pr. All in all, 
increase of B 2 , decrease v (~ Vi), and increase of I 



make the right-hand side of Equation ( 11 ) roughly 
constant, which results in a roughly linear growth 
rate of B 2 in our simulations. 

3. Third, after the slower growth stage, the saturation 
stage begins. In Run 512-B mQa; 0.0TPr, it begins at 
t~17. 

During the exponential growth stage, we can write 
B 2 ~ B 2 (t — 0)exp(f/r), where r is approximately a 
constant times the large-eddy turnover tim^J There- 
fore, the growth time of magnetic field depends on the 
value of B max . Indeed, the figure shows that the run with 
B m ax = 1 reaches the saturation stage more quickly than 
smaller B max cases. Run 512-B max 0.01-Pr (thick dashed 
curves) reaches the saturation stage later than Run 256- 
B m(M 0.0TPr (thin dashed curves), because the size of 
the initially magnetized region normalized by the system 
size is smaller in the former, which makes B 2 (t = 0) 
smaller in the former. If the seed magnetic field is ex- 
tremely weak, it is possible that the growth of magnetic 
field is slower in high magnetic Prandtl number turbu- 
lence than in the unit Prandtl number cases with the 
same seed magnetic field. 

As in the unit Prandtl number cases, the seed magnetic 
fields fills the computational box very fast. According to 
Figure [7](d) , it takes less than~ 4 large-eddy turnover 
times for the seed fields to fill the whole box. Roughly 
speaking, the standard deviation of magnetic field dis- 
tribution, ct, grows exponentially when B max is weak. 
Consider a magnetized region of size ct with a < L. 
Then, since stretching happens near the outer-scale of 
turbulence, the timescale for a substantial change of ct 
is roughly the large-eddy turnover time L/v. Therefore, 
wc have the relation 



Act ct 

At L/v 



(12) 



which implies exponential growth of ct when a < L. The 
inset of Figure [n d) , in which the vertical axis is drawn 
in logarithmic scale, supports exponential growth of ct. 

3.3. Decaying turbulence 

If the ICM turbulence is driven by temporarily inter- 
mittent events, such as mergers, the ICM will undergo 
disruption and relaxation repeatedly. In this subsection, 

3 In our simulations, r might be the eddy turnover time at 
k ~10 (see the previous footnote) 



we study growth of a localized seed magnetic field in tur- 
bulence with episodic driving. 

To see if a single merger can disperse a localized seed 
magnetic field, we first perform a decaying turbulence 
simulation (Run 256-B ma:r 0.01-f0). At t=0, we have a 
fully developed turbulence and a localized seed magnetic 
field with B ma:E =0.01. The shape of the seed magnetic 
field is given in Equation Q . Then we let the turbulence 
decay without forcing it. The solid curve in the left panel 
of Figure M shows that the standard deviation ct of the 
magnetic field distribution gradually rises and converges 
to the value of 0.5 within ~10 initial large-eddy turnover 
times (L/vt—o). Therefore, we can conclude that even a 
decaying turbulence can make a localized seed field fill 
the whole system if the driving scale of turbulence is 
comparable to the system size. The inset in Figure [H] 
shows that v 2 immediately decreases as turbulence de- 
cays, while B 2 rises for ~ 10 large-eddy turnover times 
and then levels off. As a result, the ratio of magnetic to 
kinetic energy densities increases. 

The results of decaying turbulence simulation imply 
that a single merger (with magnetic field injection at 
the end of the merger) would be enough to magnetize 
the whole ICM. Consider a localized seed magnetic field 
in a cluster with driving scale of ^400 kpc and initial 
velocity dispersion of ~ 400 km/sec. Then the large- 
eddy turnover time is ~ 10 9 years. Our result in Figure 
[8] suggests that the decaying turbulence can magnetize 
surrounding ~1 Mpc region of the cluster within the age 
of the universe. 

Of course, a faster magnetization is possible when we 
apply additional forcing. For example, we may allow for 
short-lasting driving after we inject a localized seed mag- 
netic field, which mimics a seed magnetic field injection 
during a merger. We perform a simulation to see the 
effect of additional forcing (Run 256-B maa; 0.01-f3). The 
numerical setups are identical to those of the decaying 
turbulence simulation. However, in this case we drive 
turbulence for 3 code time units, which is slightly larger 
than one large-eddy turnover time L/v ~ 2.5 code time 
units. After t = 3 code time units, we turn off forcing 
and let the turbulence decay. We plot the result in the 
left panel of Figure [8] (see the dashed curve), where we 
can see that ct rises to a value very close to the saturation 
value of 0.5 within a couple of large-eddy turnover times. 
Note that the strength of turbulence is not a crucial fac- 
tor here. Whatever the strength of the turbulence is, 
a decaying turbulence can virtually seed the whole sys- 
tem within a couple of large-eddy turnover times. But, 
if the strength of turbulence is very weak, the large-eddy 
turnover time can be very long. The inset shows that 
magnetic energy density rises for ~ 10 code time units 
and then levels off, while kinetic energy density begins 
to decrease right after we turn off forcing at t = 3 (in 
code units) . The magnetic energy density in this simula- 
tion (dashed curve) is slightly higher than that of purely 
decaying turbulence (solid curve). From this simulation 
and the purely decaying turbulence simulation, we can 
conclude that even a single merger is marginally enough 
to seed the whole cluster. However, the resulting mag- 
netic field in this case may not be so strong (see values 
of B 2 in the inset). 

The right panel of Figure [8] shows time evolution of 




Fig. 7. — High magnetic Prandtl number turbulence (y 2> rj). (a) Time evolution of B 2 and v 2 . In all simulations, the Reynolds number 
is ~ 100 and the magnetic Reynolds number is very large because we use a hyper-diffusivity. Except in case of Bo = 0.0001, where we use 
a uniform seed magnetic field, we use a localized seed magnetic field. Time is given in units of the large-eddy turnover time defined by 
L/v, where L (~ 2.5) is the outer-scale of turbulence and v (~ V0.9) is the r.m.s. velocity before the saturation stage. We observe fast 
growth of magnetic field. Results are from Runs 512-B maI 0.01-Pr, 256-B m(l:E 0.01-Pr, 256-B maa; l-Pr, and REF256-Pr. (b) Time evolution 
of velocity and magnetic spectra. Note that the velocity spectrum does not change much during the exponential growth stage. From Run 
512-B maI 0.01-Pr. (c) Same as (a), but the vertical axis is drawn in logarithmic scale. Note that all simulations with weak B max have 
similar growth rates during the exponential stage. This is because the stretching scales are virtually same in the simulations. In Run 
512-B maa; 0.01-Pr, the exponential growth stage ends at t ~12. (d) Standard deviation of magnetic field distribution. The localized seed 
magnetic field fills the whole box within ~3 to 4 large-eddy turnover times. 



magnetic spectrum in Run 256-B maa; 0.01-f3, which cor- 
responds to the dashed curves in the left paneQ When 
t < 7.5L/vt=o, magnetic energy spectrum peaks near the 
dissipation-scale wavenumber, which implies stretching is 
most active there, and the magnetic spectrum moves up- 
ward without changing its shape much. This behavior 
of magnetic spectrum is very similar to that of a driven 
turbulence. However, there are also differences^] First, 
the dissipation scale in the decaying turbulence gradually 
increases as turbulence decays. In incompressible hydro- 

4 The time evolution of magnetic spectra in Run 256-B maI 0.01- 
fO is qualitatively similar. The reason we show spectra of Run 256- 
B maa; 0.01-f3 is that they evolve more quickly due to the additional 
forcing. 

5 Since magnetic field fills most of the computational box rea- 
sonably fast (see the left panel) and evolution of the system after 
magnetic field fills the box is very similar to that of a weak uniform 
seed field case, the discussion below should be also applicable to 
the case with a weak uniform seed field. 



dynamic turbulence, the dissipation scale is ~ (Lv/v) 3 ^ 4 
times smaller than L. Therefore, as v goes down, the 
dissipation scale goes up. The eddy turnover time at the 
dissipation scale increases as v decreases and the dissipa- 
tion scale increases. Second, as the eddy turnover time 
at the stretching scale (i.e. at the dissipation scale in this 
case) changes, the growth rate of magnetic energy den- 
sity deviates from an exponential one, which is clearly 
observed in the left panel of Figure M At t ~ 7.5L/vt=o, 
the magnetic spectrum 'touches' the velocity spectrum 
at the dissipation scale. 

After t ~ 7.5L/v t= Q, the evolution of magnetic spec- 
trum is somewhat complicated. When we compare mag- 
netic spectra at t=7.5 (dashed curves) and 24 (long- 
dashed curve), we find that magnetic spectrum for large 
k values goes down and that for small k values goes up 
after t ~ 7.5L/vt=o- The spectrum for large k values de- 
creases as the velocity spectrum decreases due to decay 
of turbulence. On the other hand, the spectrum for small 
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Fig. 8. — Purely decaying turbulence (Run 256-B maI 0.01-fO; solid curves) and turbulence with short-lasting driving (Run 256-B maI 0.01- 
f3; dashed curves). In the former, we let the turbulence decay at t=0. In the latter, we drive turbulence for t < 3 and we stop driving 
after t = 3. (Left) Time evolution of energy densities (actually, B 2 and v 2 ) and the standard deviation of magnetic field distribution. 
Note that magnetic field spreads and amplifies even in purely decaying turbulence. (Right) Time evolution of energy spectra (from Run 
256-B maa; 0.01-f3). Time is given in units of L/v t —o (~ 2.5 in code units), where Vt=o (~ 1) is the velocity at t=0. 




Fig. 9. 



20 30 40 

Tim e [ = t code /(L/v t=0 )] 

Turbulence with episodic forcing (Run 256-B 



c 0.01-f e pi). We repeat forcing (3 code time units) and no-forcing (7 code time 



units) every 10 code time units. (Left) Time evolution of energy densities (actually, B 2 and v 2 ) and the standard deviation of magnetic 
field distribution. We observe efficient growth and diffusion of magnetic field. (Right) Time evolution of energy spectra. 

other simulations, we have a fully developed turbulence 
and a localized seed magnetic field at t=0. In the sim- 
ulation, we drive turbulence for 3 code time units and 
let it decay for 7 code time units. Then we repeat the 
same thing periodically. From the figure, it is evident 
that a 'saturation stage' is reached at t ~ 50 in code 
units. After t ~ 50 code time units, v 2 and B 2 do not 
show a systematic decrease or increase and just fluctu- 
ate around their average values. The inset shows that 
the standard deviation at t=10 (in code units) is slightly 
smaller than the saturation value of 0.5. But the second 
episode of driving starting at t=10 (in code units) quickly 
boosts it to the saturation value. The overall behavior of 
B 2 and v 2 is very similar to that in continuous-driving 
cases. Note that, during the decay stages after t ~ 50 
(in code units), magnetic energy density sometimes ex- 



k values increases because of stretching at large scales. 
As a result, the ratio B 2 /v 2 increases and the peak of the 
magnetic spectrum moves to larger scales. This behavior 
of Eb(k) will ultimately stop when the ratio B 2 jv 2 be- 
comes large enough so that stretching at the outer scale 
of turbulence becomes suppressed. Note that, since mag- 
netic spectrum for small k values increases and that for 
large k values decreases, the total magnetic energy den- 
sity changes more slowly than the kinetic energy density 
(see the dashed curves for t >7.5 in the left panel). 

3.4. Turbulence with episodic driving 

In order to have strong magnetization, we need either a 
continuous or a frequent episodic forcing. The left panel 
of Figure [9] shows that, if there are repeated mergers, 
we can have a relatively strong magnetic field. As in 
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ceeds kinetic energy density, which is not observed in 
continuously driven turbulence (with v — rj). Therefore, 
it may be used to distinguish the driving mechanisms in 
the ICM. 

The right panel of Figure [9] shows time evolution of 
energy spectra. During the growth stages (i.e. t <50 
in code units), the overall behavior of magnetic spec- 
trum is very similar to that in continuously driven cases. 
For example, during the rapid growth stage (t <15 in 
code units), magnetic spectrum peaks near the dissipa- 
tion scale and moves upward without changing the shape 
much (see solid and dotted curves that represent Eb(k)). 
This behavior of magnetic spectrum is very similar to 
that during the exponential growth stage in a contin- 
uously driven turbulence. The peak of magnetic spec- 
trum moves to smaller wavenumbers during the slower 
growth stage (15< t <50), which is similar to the be- 
havior of magnetic spectrum during the 'linear' growth 
stage in a continuously driven turbulence. However, after 
t >50, magnetic and velocity spectra during the driving 
episodes are clearly different from those during the decay 
(i.e. no-forcing) phases. Let us compare spectra at t =90 
(dashed curves) and 93 (long-dashed curves). Since an 
episode of driving starts at t=90 and ends in t=93 in 
code units, the dashed curves are spectra during the de- 
cay phase and the long-dashed curves are spectra at the 
end of the driving episode. Velocity spectrum during the 
decay phase (dashed curve) is clearly lower than that 
during the driving episode (long-dashed curve). Mag- 
netic spectrum during the decay phase is also lower than 
that during the driving episode. 

4. DISCUSSION 

In this paper, we have found that growth of a local- 
ized seed magnetic field is as fast as that of a uniform 
magnetic field. This result implies that magnetic field 
ejected from astrophysical bodies can be a viable origin 
of magnetic field in the large-scale structure of the uni- 
verse. 

Assuming a uniform seed magnetic field (and tur- 
bulence generated by cosmological shocks), Ryu et al. 
(2008) showed that turbulence in the ICM has almost 
reached the saturation stage. Our current simulations 
imply that we may arrive at a similar conclusion also in 
the case of a localized seed magnetic field. Therefore, we 
expect a strong magnetization in clusters regardless of 
the shape of the seed field, which makes the findings in 
Ryu et al. (2008) more reliable in clusters. 

Since both a uniform and a localized seed fields give 
similar distributions and strengths of magnetic fields in 
clusters, it may be difficult to tell which seed field is 
more important. One may think that it would be more 
advantageous to look at magnetic field distributions in 
filaments. Indeed, we expect that large-eddy motions 
are slower in filaments and magnetic fields ejected from 
astrophysical bodies may not have enough time to fill the 
whole system (see also discussions in Bruggen et al 2005) . 
Therefore, magnetic field distribution in filaments may 
be spatially inhomogeneous when seed magnetic fields are 
ejected from galaxies. Indeed, if t/(L/v) < 2, we expect 
to see spatially inhomogeneous distribution of magnetic 
fields in filaments (see Figures [5] and [6] see also the right 
panel of Figure 0. Even if tf{L/ v) is larger, driving 
scale of turbulence can also affect spatial inhomogeneity. 



If turbulence in filaments is driven on scales much smaller 
than the system size of a few Mpc, magnetic fields may 
not have enough time to diffuse over un-correlated eddies 
and, as a result, we may have spatially inhomogeneous 
magnetic field distributions. 

On the other hand, if turbulence in a filament is driven 
by cosmological shocks, where driving scale of turbulence 
is comparable to the system size, we will have more or less 
homogeneous magnetic field distribution. Cho & Ryu 
(2009) argued that the driving scale of turbulence is com- 
parable to the size of the system and that t/(L/v) ~ 4 
in a typical filament]^] Figure M suggests that 4 large- 
eddy turnover times are more than enough to magnetize 
the whole outer-scale eddy and also adjacent outer-scale 
eddies, the total size of which is larger than a few Mpc. 
Therefore, unless the astrophysical sources that provide 
seed magnetic fields are very sparse, magnetic fields fill 
substantial volume fraction in filaments. 

Our results show that magnetic diffusion is fast for 
B ma a; _ 1- Note that, when T5 max — 1, local magnetic 
energy density is comparable to kinetic energy density 
and magnetic backreaction is not negligible in the local 
region. Nevertheless, our simulations show overall mag- 
netic diffusion by turbulence is not suppressed by the 
backreaction. In this paper, we do not consider stronger 
initial seed magnetic fields because we can guess what 
will happen. In case of B ma2 , > 1, we expect that mag- 
netic pressure first makes magnetic field spread outward, 
which should be faster than turbulence diffusion (see Col- 
bert et al. 1996; Xu et al. 2010, 2011). Then, after local 
magnetic energy density drops below the kinetic energy 
density, turbulence becomes the major agent that makes 
magnetic field spread out, which we know from our cur- 
rent simulations is fast. Therefore we expect fast diffu- 
sion of the seed magnetic field even cases of B max > 1. 

We do not assume any specific origin of turbulence 
in this paper: unless the driving scale is much smaller 
than the system size, it does not matter much. Con- 
sider a cluster of galaxies. If turbulence is initiated by 
cosmological shocks (Ryu et al. 2003; Pfrommer 2006) or 
mergers (De Young 1992; Tribble 1993; Norman & Bryan 
1999; Roettiger, Burns, & Stone 1999; Ricker & Sarazin 
2001), we expect that the outer scale of turbulence is 
comparable to the size of the whole system. In fact, the 
outer scale of turbulence observed in some simulations 
during the formation of galaxy clusters is up to ~ sev- 
eral 100 kpc (Norman & Bryan 1999; Ricker & Sarazin 
2001), which is a few times smaller than the cluster size 
of ~ Mpc. In our current simulations, the driving scale is 
about 2.5 times smaller than the simulation box. There- 
fore, our current simulations mimic turbulence generated 
by cosmological shocks or mergers. However, it is still 
possible that turbulence in the ICM is driven by infall 
of small structures (Takizawa 2005), AGN jets (see, for 
example, Scannapieco & Bruggen 2008), or galaxy wakes 
(Roland 1981; Bregman & David 1989; Kim 2007). In 
this case, the driving scale of turbulence can be much 
smaller than the size of the whole system. If we consider a 

6 Cho & Ryu (2009) used a different definition of the large-eddy 
turnover time. The large-eddy turnover time t et idy m Cho & Ryu 
(2009) is roughly 2.5-3 times smaller than our large-eddy turnover 
time. Cho & Ryu (2009) found that t/t e ^dy ~ 10 f° r a typical 
filament, which is equivalent to ~4 large-eddy turnover times in 
our definition here. 
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localized seed magnetic field inside an energy-containing 
eddy, we know that the seed field can spread and fill the 
outer-scale eddy fast. If the driving scale of turbulence 
is small, the subsequent stage of magnetic field diffusion 
may take a long time to fill the whole system because dif- 
fusion of magnetic field over un-correlated eddies may be 
slower compared with that inside an energy-containing 
eddy. We will address this issue elsewhere (Cho 2012, in 
preparation) . 

In this paper, we have considered homogeneous turbu- 
lence. However, turbulence in the ICM is not homoge- 
neous. The ICM density varies from the core region to 
the cluster outskirts. Then how can it affect our results? 
The value of the ICM density length-scale may be an 
important factoij^ If the length-scale is larger than the 
outer scale of turbulence, our results in this paper will 
not be affected much. On the other hand, if the opposite 
case is true, then it depends on the driving mechanism 
in the following way. First, if turbulence is initiated by 
a violent event /events and if the outer scale of turbu- 
lence does not change for a couple of large-eddy turnover 
times, then our results may not be affected much. If the 
outer scale of turbulence does not change, a couple of 
large-eddy turnover times would be enough for magnetic 
field to fill the most of the ICM. Second, if turbulence is 
initiated less violently, then it may be possible that the 
effective outer scale of turbulence reduces to a value com- 
parable to the density length-scale. In this case, effective 
reduction of the outer scale can make the resulting tur- 
bulence similar to that with small-scale driving. Further 
numerical study will be necessary to test this possibility. 

5. SUMMARY 

In this paper, we have assumed that driving scale of 
turbulence is comparable to the size of the system, which 
can be justified if we consider ICM turbulence driven 
by mergers or cosmological shocks. We have found the 
following results: 

1. The growth of a localized seed magnetic field is as 



fast as that of a uniform seed magnetic field. 

2. Turbulence diffusion of magnetic field is fast. When 
we insert a localized seed magnetic field in a tur- 
bulent medium, turbulent motions make it disperse 
and fill the whole system very fast, within ^3 large- 
eddy turnover times. After magnetic field fills the 
whole system, the time evolution should be very 
similar to the case of a uniform seed magnetic field. 
This way, we have a fast growth of the seed mag- 
netic field. Our result implies that even in fila- 
ments the volume filling factor of magnetic fields 
is of order unity if turbulence is driven on scales 
comparable to the size of the systems. 

3. Growth and turbulence diffusion of a localized seed 
magnetic field in a high Prandtl number turbulence 
is also fast. 

4. A localized seed magnetic field can ultimately fill 
the whole system even in a decaying turbulence. 
Therefore, we may have a wide-spread magnetic 
field in a cluster that has undergone just one ma- 
jor merger. Magnetic field can be injected at any 
time during the merger: our decaying turbulence 
simulation suggests that magnetic field can be in- 
jected even at the end of the merger. Although the 
strength of the resulting magnetic field is very weak 
in decaying turbulence, the wide-spread magnetic 
field can be used as a seed field in the next episode 
of merger. 

5. Growth of a localized seed magnetic field in case of 
episodic driving is also fast. 



This research was supported by National R&D Pro- 
gram through the National Research Foundation of Ko- 
rea (NRF) funded by the Ministry of Education, Sci- 
ence and Technology (No. 2011-0018751). We thank the 
anonymous referee for constructive comments. 



REFERENCES 



Arieli, Y., Rephaeli, Y. & Norman, M. L. 2011, ApJ, 738, 15 
Banerjee, R., & Jedamzik, K. 2003, Phys. Rev. Lett., 91, 251301 
Batchelor, G. 1950, Proc. R. Soc. London A, 201, 405 
Bertone, S., Vogt, C, & Enfilin, T. 2006, MNRAS, 370, 319 
Brandenburg, A. & Subramanian, K. 2005, Phys. Reports, 417, 1 
Bregman, J. N., & David, L. P. 1989, ApJ, 341, 49 
Briiggen, M., Ruszkowski, M., Simionescu, A., Hoeft, M., & Dalla 

Vecchia, C. 2005, ApJ, 631, L21 
Childress, S., & Gilbert, A. 1995, Stretch, Twist, Fold: The Fast 

Dynamo (Berlin: Springer) 
Cho, J., Lazarian, A., Honein, A., et al. 2003, ApJ, 589, L77 
Cho, J., & Vishniac, E. T. 2000, ApJ, 538, 217 
Cho, J., Vishniac, E. T., Beresnyak, A., Lazarian, A., & Ryu, D. 

2009, ApJ, 693, 1449 
Cho, J. & Ryu, D. 2009, ApJL, 705, 90 

Colbert, E. J. M., Baum, S. A., Gallimore, J. F., O'Dea, C. P., & 

Christensen, J. A. 1996, ApJ, 467, 551 
De Young, D. 1992, ApJ, 386, 464 

Dolag, K., Bartelmann, M., & Lesch, H. 1999, A&A, 348, 351 
Dolag, K., Bartelmann, M., & Lesch, H. 2002, A&A, 387, 383 

7 In addition, magnetic buoyancy can also be an important 
factor. If a localized seed magnetic field is injected into the central 
region of the ICM, magnetic buoyancy will enhance the rate of 
magnetic diffusion. 



Donnert, J., Dolag, K., Lesch, H., Miilller, E. 2009, MNRAS, 392, 
1008 

Dubois, Y., & Teyssier, R. 2008, A&A, 482, L13 
Ensslin, T. A., Biermann, P. L., Kronberg, P. P., & Wu X. 1997, 
ApJ, 477, 560 

Gnedin, N. Y., Ferrara, A., & Zweibel, E. G. 2000, ApJ, 539, 505 
Kazantsev, A. P. 1968, Soviet Phys.-JETP Lett., 26, 1031 
Kim, W.-T. 2007, ApJL, 667, 5 
Kronberg, P. P., 1994, Rep. Prog. Phys., 57, 325 
Kronberg, P. P., Lesch, H., & Hopp, U. 1999, ApJ, 511, 56 
Kulsrud, R. & Anderson, S. 1992, ApJ, 396, 606 
Kulsrud, R. M., Cen, R., Ostriker, J. P., & Ryu, D. 1997, ApJ, 
480, 481 

Kulsrud, R. M., & Zweibel, E. G. 2008, Reports on Progress in 

Physics, 71, 046901 
Lesieur, M. 1990, in Turbulence in Fluids: Stochastic and 

Numerical Modelling (2d rev. ed.; Dordrecht: Kluwer), chap. 8 
Murgia, M., Govoni, F., Feretti, L., Giovannini, G., Dallacasa, D., 

Fanti, R., Taylor, G. B., & Dolag, K. 2004, A&A, 424, 429 
Norman, M. L., & Bryan, G. L. 1999, Lecture Notes in Physics 

Vol. 530: The Radio Galaxy Messier 87, 530, 106 
Pfrommer, C, Springel, V., EnBlin, T. A., & Jubelgas, M. 2006, 

MNRAS, 367, 113 
Rees, M. 1987, QJRAS, 28, 197 

Rephaeli, Y. 1988, Comments on Astrophysics, 12, 265 



Growth of Magnetic Field 



11 



TABLE 1 
Simulations. 



Run 


Resolution 


$Tnax 3<t t — 


B h 


Pr,m c 


forcing 


256-B max 0.001 


256 3 


0.001 


0. 


1 


continuous 


256-B ma:c 0.01 


256 3 


0.01 


0. 


1 


continuous 




256 3 


0.1 


o 


I 


r'AVltlTlllAIIC 
LWIILIIILIW U.S 


256-B ma3; 1 


256 3 


1 


0. 


1 


continuous 


512-B maa; 0.01 


512 3 


0.01 


0. 


1 


continuous 


256-B ma:E 0.01-Pr 


256 3 


0.01 


0. 


high 


continuous 


256-B ma:E l-Pr 


256 3 


1 


0. 


high 


continuous 


512-B max 0.Ql-Pr 


512 3 


0.01 


0. 


high 


continuous 


256-B maa; 0.01-fO 


256 3 


0.01 





1 


no forcing 


256-B maa; 0.0l-f3 


256 3 


0.01 





1 


t < 3 in code units 


256-B maa; 0.01-f (!pi 


256 3 


0.01 





1 


episodic 


REF256 


256 3 




10~ 4 


1 


continuous 


REF256-Pr 


256 3 




io- 4 


high 


continuous 
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